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ABSTRACT 


The  atmospheric  range  error  in  an  electromagnetic 
signal  arriving  vertically  at  the  earth  is  measured  by  the 
height  integral  of  the  refractivity  N  through  the  atmosphere. 
(N  =  10®  (n-1)  where  n  is  the  index  of  refraction. )  N  in 
the  non-ionized  atmosphere  is  the  sum  of  "dry"  and  "wet" 
components,  which  are  considered  separately.  Only  the 
dry  part  is  important  at  optical  wavelengths,and  the  present 
discussion  deals  principally  with  this  dry  part.  Data  from 
several  thousand  meteorological  balloons  were  studied  to 
develop  a  method  for  predicting  the  height  integrals  of 
Ndry*  ft  is  found  theoretically  and  from  observed  data 
that, although  N<jry  depends  on  both  temperature  and 
pressure,  its  height  integral  is  a  linear  function  of  surface 
pressure  only.  The  slope  is  nearly  the  same,  worldwide; 
there  is,  however,  a  very  small  but  detectable  latitude 
variation,  corresponding  to  the  latitude  variation  of  g.  It 
is  shown  that  the  vertical  range  error  (dry  component)  at 
a  given  latitude  can  be  predicted  within  2  mm  rms  from 
surface  pressure  alone.  Irregularities  in  the  temperature 
profile  do  not  affect  the  prediction  in  the  zenith  direction. 

In  view  of  these  results,  the  equivalent  height  for  a  quartic 
model  must  vary  as  surface  temperature;  it  is 
approximately  40  km  at  0°C  and  becomes  zero  near  0°K  at 
all  locations.  Possible  small  differences  between 
parameters  for  different  stations  are  being  examined 
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1 .  INTRODUCTION 


The  electromagnetically  measured  range  from  the 
earth  to  an  object  outside  the  atmosphere,  e.  g. ,  the  moon 
or  an  artificial  satellite,  is  too  large  by  more  than  2  meters 
if  the  object  is  directly  overhead,  and  by  about  twice  that 
amount  if  the  zenith  angle  of  the  line  of  sight  is  60°.  The 
precise  amount  of  the  atmospheric  effect  on  range  is  given 
by  the  integral  of  (n-1)  along  the  signal  path,  where  n  is 
the  index  of  refraction.  For  convenience,  the  refractivity 
N  is  defined  as  10b  (n-1).  Since  N  depends  on  atmos¬ 
pheric  conditions,  it  varies  along  the  signal  path. 

In  the  optical  region,  N  for  air  is  also  a  function 
of  wavelength,  and  range  measurements  at  two  optical 
wavelengths  could  be  used  to  remove  refraction  effects 
from  the  data.  The  instrumentation  requirements  for  this 
are,  however,  very  severe.  In  the  radio  region,  the  re¬ 
fractivity  of  the  uncharged  atmosphere  (troposphere  and 
stratosphere)  is  independent  of  signal  frequency  up  to 
15  000  MHz. 

A  theoretical  refraction  correction,  if  adequate, 
would  be  simpler  than  multifrequency  measurements.  A 
model  to  provide  this  correction  was  developed  for  radio 
use  (initially  for  use  with  satellite  doppler  data),  and  the 
so-called  "dry"  part  of  it  can  also  be  applied  in  the  optical 
region  with  little  change. 
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2.  THEORY 


In  order  to  predict  the  atmospheric  range  error  in 
an  electromagnetic  signal  at  any  elevation  angle,  it  is 
necessary  to  know  the  range  error  for  a  signal  in  the  zenith 
direction,  i.  e. ,  the  height  integral  10“6jNdh  in  the  exist¬ 
ing  state  of  the  atmosphere.  The  error  at  a  small  zenith 
angle  can  be  obtained  directly  from  this,  being  little  affected 
by  the  height  profile  of  N.  For  very  large  zenith  angles, 
however,  it  is  necessary  to  know  the  N  profile  also. 


ZENITH  INTEGRAL  OF  N 

The  refractivity  of  air  at  radio  frequencies  up  to 
15  000  MHz  can  be  expressed  as  the  sum  of  two  components, 
"dry"  and  "wet,  "  which  will  be  designated  by  the  subscripts 
d  and  w,  respectively.  Their  values  at  a  given  point  of  the 
atmosphere  can  be  computed  from  meteorological  data.  The 
equations  given  by  Smith  and  Weintraub  (Ref.  1)  are: 

77.  6  P  ) 


3.  73  x  10  e 


in  which  T  is  the  temperature  in  degrees  Kelvin,  P  is 
the  total  pressure  in  millibars,  and  e  is  the  partial  pres¬ 
sure  of  water  vapor,  also  in  millibars. 

At  optical  frequencies,  the  effect  of  water  vapor  is  so 
much  decreased  as  to  be  practically  negligible.  The  above 
expression  for  N^,  however,  is  in  close  agreement  with 
values  of  N  given  by  the  Edlen  expression  (Ref.  2)  for  the 
phase  refractivity  of  dry  air  in  the  infrared  (wavelengths  of 
a  few  microns).  The  following  discussion  of  the  height 
integral  of  N^  therefore  pertains  either  to  the  dry  component 
of  radio  range  error  or  to  the  entire  optical  range  error  in 
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the  near  infrared;  but  in  the  optical  case,  a  dispersion 
correction  of  a  few  percent  will  be  needed  for  laser 
wavelengths  in  the  visible  (Ref.  3). 

The  height  integral  of  in  Eq.  (1)  can  be  deduced 
from  the  gas  laws  without  knowledge  of  the  shape  of  the  Nd 
profile.  If  V  is  the  volume  and  R  is  the  gas  constant 
for  unit  mass  of  dry  air,  if  we  assume  that  air  is  a  perfect 
gas,  we  have,  from  the  gas  laws, 

PV  =  RT 

|  =  f  =  R P  »  (2) 


where  p  is  the  density.  The  combination  of  Eqs.  (1)  and 
(2)  yields: 


N.  =  10  x  77.  6  Rp  , 
d 


if  R  and  p  are  in  cgs  units.  Then 


/ 


N  .  dh  =  10"3  x  77 
d 


.enjpdh 


(3) 


The  pressure  at  the  surface  is: 


=  IP  g  dh 


(4) 


If  g,  the  acceleration  of  gravity,  is  assumed  con¬ 
stant  over  the  height  range  of  the  lower  atmosphere,  the 
combination  of  Eqs.  (3)  and  (4)  yields  an  atmospheric 
height  contribution  of: 


Ah 


/ 


77.  6RP 

(n-l)dh  =  10"  x  ■ 


g 


(5) 
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if  all  quantities  are  in  cgs  units.  Thus  theoretically, in  a 
dry  atmosphere,the  zenith  range  error  produced  by  refrac¬ 
tion  is  directly  proportional  to  the  surface  pressure  and  is 
independent  of  the  way  in  which  temperature  varies  with 
height  (Refs.  4  and  5).  Observed  data  to  test  this  conclusion 
for  the  real  atmosphere  will  be  presented  below. 

The  above  theory  does  not  apply  to  the  wet  component, 
which  will  be  discussed  briefly  later. 


HEIGHT  PROFILE  OF  N 


Although  the  N  profile  is  not  needed  to  predict  the 
atmospheric  effect  at  or  near  the  zenith,  this  is  not  true  at 
low  elevation  angles.  To  be  theoretically  valid,  a  profile 
model  must  account  for  a  zenith  integral  of  N  which  is 
consistent  with  Eq.  (5).  In  developing  a  model  it  is  assumed 
that  each  component  of  N  (dry  or  wet)  is  a  function  of  its 
surface  value  and  of  height  above  the  earth,  but  not  of  hori¬ 
zontal  position  or  time  during  a  short  series  of  observations. 
Thus  diurnal  variations  and  weather  fronts  are  neglected 
for  the  time  being.  The  dry  and  wet  components  are 
different  functions  of  height. 


In  a  dry  atmosphere,  at  heights  small  enough  so  that 
the  acceleration  of  gravity  g  can  be  considered  constant, 
the  pressure  (Ref.  6)  is  given  by: 


P 


P  e 
s 


.6.  fhSh 

R  Jo  T 


(6) 


where  Pg  is  the  pressure  at  the  surface  of  the  earth,  R 
is  the  gas  constant, and  h  is  height  above  the  surface. 


If  the  temperature  T  does  not  vary  with  height  in 
the  atmosphere,  this  becomes: 

.  JSlL 
RT_ 


P  = 


P  e 
s 


(7) 


-  5  - 


TNI  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
Silver  Skkino  Maryland 


Combining  Eq.  (7)  with  Eq.  (1)  yields  an  exponential 
N  profile  for  a  dry  isothermal  atmosphere: 

RT 

N  =  N  e  8  .  (8) 

s 


Exponential  models  for  the  total  N  and  for  its 
separate  components  have  been  extensively  used  (Refs.  7, 
8,  and  9). 


In  the  earth's  atmosphere,  however,  the  temperature 
on  the  average  decreases  with  height  at  a  fairly  constant 
lapse  rate  in  the  troposphere,  is  moderately  constant  in  the 
tropopause  region,  and  above  that  increases  slowly  with 
height  in  the  stratosphere.  Considering  the  troposphere 
first,  let  us  assume  that  the  lapse  rate  a  is  constant 


(° 5  - 1)  •  Then: 


T  = 


-  ah  . 


(9) 


Substituting  Eq.  (9)  into  Eq.  (6)  and  integrating, 
we  get  (Ref.  6) : 


P  = 


(10) 


Thus,  the  pressure  when  a  ^  0  is  not  an  exponential 
function  of  height.  If  P  and  T  from  Eqs.  (9)  and  (10)  are 
used  in  Eq.  (1),  the  profile  can  be  put  in  the  form 
(Ref.  4):  _  Q 

Tt  1  w 


(U) 
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where 


i» 


(12) 


The  degree  of  the  profile  is  shown  in  Fig.  1  as 
a  function  of  the  temperature  lapse  rate  a.  For  any  posi¬ 
tive  value  of  a  (negative  -aX  )  which  is  likely  to  occur  in 
the  troposphere,  4  is  positive. 


The  degree  n  would  be  zero  for  a  temperature  lapse 
rate  of  34°C/km,  an  improbably  high  value.  The  adiabatic 
lapse  rate  is  9.  8°C/km  (Ref.  6),  but  observed  tropospheric 
lapse  rates  are  generally  less  (~  7°C/km  in  warm  climates, 
less  near  the  poles),  a  and  M.  are  negative  in  the  region 
of  a  temperature  inversion,  and  negative  in  the  stratosphere 
above  the  tropopause. 

It  was  shown  earlier  (Ref.  5)  that  the  zenith  integral 
of  Nj  in  an  atmosphere  with  constant  lapse  rate  (zero, 
positive  or  negative),  as  obtained  by  integrating  Eq.  (8) 
or  Eq.  (11)  with  appropriate  limits,  is  identical  to  the  value 
of  Eq.  (5).  If  the  lapse  rate  is  constant  and  known,  it  is 
therefore  possible  to  write  an  equation  for  the  profile 
which  will  match  an  observed  profile  as  regards  both  the 
zenith  integral  and  the  profile  shape.  Since  in  practice  the 
lapse  rate  is  only  fairly  constant  in  the  troposphere  and 
changes  sign  at  the  tropopause,  any  single  mathematical 
function  will  not  match  the  actual  profile  shape  per¬ 
fectly  at  all  heights.  A  single  function  can,  however,  yield 
the  correct  zenith  integral,  regardless  of  irregularities  in 
the  profile  shape;  and  also  provide  a  usable  approximation 
to  the  profile  shape  in  the  denser  part  of  the  atmosphere. 

A  fourth-degree  model  (|i  =  4  in  Eqs.  (11)  and  (12)), 
which  was  developed  earlier  (Ref.  4),  corresponds,  if 
g  =  980  cm/s  ,  to  a  temperature  lapse  rate  of  6.  828°C/km, 
and  does  on  the  average  match  observed  profiles  well  to  a 
considerable  height,  in  regions  of  the  earth  where  this 
lapse  rate  is  a  realistic  value  in  the  troposphere  (Fig.  2). 
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Fig.  1  DEGREE  OF  THEORETICAL  Nd  PROFILE 
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Fifl.  2  PROFILE  OF  TROPOSPHERIC  REFRACTIVITY  N 
WEATHER  SHIP  E,  35°|u  48°  W,  JULY  1967 
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The  quartic  model  provides  an  integrable  expression 
for  the  range  correction  JTn  dp  along  the  slant  range  vector 
p  at  any  elevation  angle,  and  this  can  be  used  as  a  range 
correction  at  elevation  angles  for  which  the  bending  of  the 
signal  path  is  negligible.  An  algorithm  for  simplifying  the 
computation  of  the  correction  has  been  developed  by 
Yionoul,is  (Ref.  10). 

Height  parameter^  for  the  quartic  model  have  been 
obtained  from  observed  data  that  will  be  described  below. 

Let  us  write  the  quartic  profile  expression  in  the  form: 

'  Nd  ■  . 

<Nd  =  — f  <hd-h>  ,  (13) 

hd 

whef*e  the  height  hd  replaces  Tg/a  of  Eq.  (11).  For 
Tg  =  273.  16°K  (0°C)  and  a  =  6.  828°/km.  T  /a  =  40.  006  km; 
but  hd  is  now  to  be  determined  empirically.  Heights  are 
to  be  measured  above  the  surface.  The  zenith  range  effect 
is  now: 

ii  ’ 

N. 

r  /*  R  dq  d 

'  Ah  =  10  INd  dh  =  10  •  — | - ,  (14) 

in,  which  h^  and  the  height  error  Ah  are  in  the  same  units. 
In  view  of  Eq.  (1),  this  can  be  compatible  with  Eq.  (5)  only 
if  h^  ;  is  proportional  to  Tk>  This  hypothesis  was  tested  by 
examining  balloon  meteorological  data,  as  will  be  described 
in  the  hext  section. 

Height  parameters  for  the  wet  component  Nw  for 
the  radio  case  were  also  obtained  from  the  data.  The 
theoretical  considerations  given  above  do  not  apply  to  the 
wet  term,  but  for  convenience,  it  has  been  assumed  that  a 
quartic  profile  .can  be  used  for  Nw  also.  The  work  on  the 
wet  component  will  be  mentioned  only  briefly;  it  is  less  far 
advanced,  and  also  is  not  of  interest  for  laser  applications. 
Except  at  very  low  angles,  the  dry  component  accounts  for 
at  least  90%  of  the  range  effect  even  at  radio  frequencies. 

! 
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The  wet  component,  however,  is  responsible  for  most  of  the 
variability  of  the  radio  range  effect  and  is  much  less 
predictable  than  the  dry  part.  Much  more  work  is  needed 
on  the  wet  component. 
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3.  DATA  AND  COMPUTATIONS 


Meteorological  balloon  data  in  14  one->ear  sets,  from 
ten  different  geographical  locations,  have  been  studied. 

There  were  two  balloons  per  day  in  each  case:  so  far,  a 
total  of  approximately  10  000  balloons.  For  each  balloon 
flight,  data  were  obtained  at  the  standard  pressure  levels 
and  also  at  additional  "significant  levels,  "  i.  e. ,  generally 
50  or  60  observed  points  per  balloon  ascent.  The  data  at 
each  point  included  height,  pressure,  temperature,  and 
relative  humidity,  thus  providing  the  necessary  information 
for  computing  the  dry  and  wet  components  of  N  from  Eq. 

(1)  (for  radio  frequencies).  Each  balloon  yielded  one  observed 
profile  of  N.  The  data  were  on  magnetic  tape  and  the  com¬ 
putations  were  done  with  an  IBM  360/91  computer. 

The  zenith  integrals  of  N^  and  Nw  were  obtained 
from  these  N  profiles  by  numerical  integration.  Occasional 
flights  in  which  any  sort  of  equipment  failure  occurred  were 
deleted  from  the  data  collection.  The  data  used  included 
only  flights  which  gave  complete  dry  data  to  at  least  thfe 
30 -mb  level  24  km)  and  complete  wet  data  to  at  least  the 
500-mb  level  (5.  5  to  6  km).  In  computing  the  numerical 
integrals,  a  small  correction  was  added  at  the  top  so  as  to 
approximate  the  true  total  integral.  For  the  dry  component, 
this  correction  was  based  on  the  topmost  observed  pressure 
and  Eq.  (5)  (a  slight  improvement  over  the  correction  used 
in  Ref.  5). 


Two  one -year  sets  of  the  zenith  integrals  obtained 
from  these  balloon  data  are  shown  in  Figs.  3  and  4:  one 
for  Washington,  D.  C.  (85  meters  above  sea  level),  and  one 
for  Albuquerque,  New  Mexico  (1620  meters  above  sea  level). 
The  three  sections  of  these  figures  show,  respectively,  the 
total  integral  and  the  separate  integrals  of  the  wet  and  the 
dry  components,  as  functions  of  time  during  the  year.  Each 
point  represents  data  from  a  single  balloon  flight. 

It  is  noticeable  that  the  dry  component  integral  shows 
little  seasonal  variation,  although  some  weather  effects  can 

Preceding  page  blank 
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be  seen.  The  wet  component  shows  large  seasonal  and  also 
weather  effects. 

The  wet  integral  J*NW  dh  is  less  than  dh  by 

an  order  of  magnitude  or  more,  yet  produces  most  of  the 
variation  in  the  total  integral,  at  radio  frequencies.  At 
laser  wavelengths,  the  wet  contribution  is  so  much  reduced 
as  to  be  practically  negligible. 

The  average  value  of  dh  through  the  atmosphere, 
for  the  year,  was  approximately  2.  3  meters  at  Washington, 
while  that  at  Albuquerque  was  40  cm  less,  or  1.  9  meters, 
because  of  the  different  heights  of  the  two  stations  above  sea 
level.  The  average  radio  value  of  Jn\v  dh  for  the  year  at 
Albuquerque  was  approximately  60"  of  that  at  Washington. 
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4.  RESULTS 


J Nd  dh  VERSUS  SURFACE  PRESSURE 

Figure  5  shows  samples  of  the  height  integral  of  Nd 
as  a  function  of  surface  pressure  for  several  samples  of 
data  from  different  locations.  Equation  (5)  predicts  a  linear 
relation.  The  plotted  points  are  observed  values  of  the 
integral  from  the  balloon  data.  It  was  assumed  (cf.  Eq.  (5)) 
that: 


/N,  dh  =  kPfl  ,  (15) 

and  the  slope  k  was  obtained  from  each  one-year  set  of 
observed  data  by  a  least  squares  procedure.  The  results 
are  listed  in  Table  1,  and  the  lines  given  by  Eq.  (15)  with 
these  empirical  slopes  k  are  shown  respectively  for 
the  samples  of  Fig.  5.  For  each  one-year  data  set,  the  rms 
error  in  predicting  dh  from  the  empirical  value  of  k 
was  between  1  and  2  mm.  The  theoretical  line  of  Eq.  (5) 
was  also  obtained  for  each  case,  using  for  each  a  mean 
value  of  g  at  mid-atmosphere  (6  km  high,  approximately 
the  500  mb  level),  at  the  latitude  of  the  station  (Ref.  11). 

The  theoretical  and  empirical  lines  are  so  nearly  coincident 
that  they  cannot  be  shown  as  separate  on  the  scale  of  the 
figure. 


The  delaying  effect  of  the  dry  atmosphere  (tropo¬ 
sphere  and  stratosphere)  on  a  ranging  signal  in  the  zenith 
direction  can  therefore  be  predicted  with  an  rms  deviation 
of  the  order  of  2  mm  of  range  (one-way  travel  of  the  signal). 

The  empirical  slopes  for  the  different  locations  are 
nearly  but  not  quite  the  same.  The  theoretical  slope  (Eq. 
(5))  is  a  function  of  the  local  value  of  g,  hence  it  varies 
slightly  with  latitude.  Figure  6  shows  this  variation.  The 
curve  represents  the  theoretical  slope  as  a  function  of 
latitude,  using  for  g  its  value  at  the  local  latitude  at  a 
height  of  6  km  above  sea  level.  The  separate  points  on  the 
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Fig.  5  (Continued) 
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Fig.  6  ZENITH  RANGE  ERROR  PER  MILLIBAR  OF  SURFACE 
PRESSURE,  DRY  COMPONENT  OF  TROPOSPHERE 
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same  graph  are  the  empirical  slopes  k  obtained  for  the 
various  data  sets;  these  also  show  a  latitude  variation  and 
their  trend  is  in  good  agreement  with  the  theoretical  variation. 
The  empirical  slope  k  for  Albuquerque  (1620  meters  high) 
is  not  markedly  different  from  the  empirical  slopes  for  the 
stations  near  sea  level. 

It  is  not  certain  that  no  biases  are  present,  but  it 
appears  that  they  must  be  small.  The  data  heights  used  in 
the  numerical  integration  are  geopotential  heights  derived 
from  the  gas  laws,  but  these  differ  very  little  from  geometric 
heights,  especially  at  middle  latitudes.  The  use  of  a  mean 
(mid-atmosphere)  value  of  g  to  get  the  theoretical  slope  of 
the  line  of  Eq.  (5)  is,  of  course,  an  approximation.  Water 
vapor  content  was  neglected  as  regards  the  value  of  the  gas 
constant  per  gram  of  air,  but  this  effect  is  vary  small.  The 
value  used  was  R  =  2.  8704  x  106  erg  g"loK_1  (Ref.  11). 

The  height  integrals  /Nd  dh  have  been  referred  to 
as  "zenith  integrals,  "  but  this  is  something  of  a  misnomer. 
Balloons  are  blown  by  the  wind  and  do  not  actually  rise 
vertically.  The  height  specified  is  the  height  co-ordinate 
of  the  balloon  position,  its  lateral  displacement  being  un¬ 
known.  The  fact  that  the  integral  of  Nj  with  the  height 
component  of  the  oblique  rise  of  the  balloon  is  so  closely 
related  to  pressure  at  the  launch  site  would  seem  to  argue 
that  the  horizontal  gradients  involved  are  not  very  large. 


HEIGHT  PARAMETERS  FOR  THE  QUARTIC  Nd  MODEL 

The  prediction  of  only  the  height  integral  of  N  may 
not  provide  enough  information  for  a  range  correction  at 
larger  zenith  angles,  where  the  temperature  profile  does 
have  a  significant  effect.  For  this,  a  model  of  the  N 
pi  ofile  may  also  be  needed. 

Since  N<j  is  proportional  to  P/T,  but  its  height 
integral  is  proportional  to  surface  pressure  alone,  it  follows 
from  Eq.  (14)  that  the  equivalent  height  hd  for  the  quartic 
Nj  profile  must  be  proportional  to  surface  temperature.  To 
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examine  this  theory,  it  was  assumed  that: 


h 


d 


+  adT 


sC  ’ 


(16) 


where  h.  is  the  "equivalent  height"  of  the  quartic 
profile  which  will  match  the  theoretical  height  integral  of 
Nd  to  the  observed  one;  h^  is  the  value  of  h^  when  the 

surface  temperature  is  0°C,  Tg^,  is  the  observed  surface 
temperature  in  °C,  and  a  is  the^temperature  coefficient 
of  the  height  increase.  The  parameters  hd  and  ad  were 
obtained  from  the  data  by  a  least-squares  procedure. 

To  do  this,  Eq.  (16)  was  combined  with  Eq.  (14)  to 
get  a  theoretical  expression  for  Jn^  dh  corresponding  to 
the  surface  Njj  at  the  time  of  a  balloon  flight;  these  theo¬ 
retical  integrals  based  on  the  observed  values  of  Ns  were 
then  equated  to  the  observed  integrals  from  balloon  data. 

The  results  are  listed  in  Table  2.  The  heights  hd  as 
listed  are  heights  above  the  station  in  each  case.  °  These 
were  found  to  be  closely  the  same  for  all  stations  regardless 
of  the  station  height  above  sea  level,  and  are  in  close  agree¬ 
ment  with  the  value  predicted  by  theory  (Ts/a,  quoted  above). 
These  heights  should  be  interpreted  as  parameters  for 
matching  observed  integrals,  not  as  indicating  any  undulation 
of  actual  pressure  levels  above  surface  undulations.  The 
results  may  contain  information  not  as  yet  deciphered, 
about  latitude,  longitude,  and  surface  height  effects.  The 
balloons  at  all  stations  were  launched  at  0  and  12  hours 
U.  T. ,  hence  at  different  values  of  local  time. 

Extrapolating  from  the  empirical  values  of  Table  2, 
the  equivalent  height  of  the  model  atmosphere  would  in  each 
case  fall  to  zero  within  a  few  degrees  of  0°K.  The  deviations 
from  absolute  zero  may  be  related  to  second-order  effects 
(e.  g. ,  gravity  anomalies,  atmospheric  tides,  etc. ).  These 
have  not  yet  been  investigated. 
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Height  Parameters  for  the  Quartic  N  .  Model  (cf.  Eq.  16) 
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!  The  integral  of,N(|  through  the  atmosphere  can  be 

predicted  at  the  zenith  or  in  a  different  direction  by  using 
the  parameters  of  Table  2  and  the  expressions  developed  in 
!  Refs.  4,  5,  and  10.  Figures  7  and  8  show  samples  of  pre¬ 
dicting  the  zenith  range  error  (dry  part)  from  the  quartic 
profile  and  the  parameters  of  Table  2.  The  rms  values  of 
i  the  deviations  are  less  than  2  mm  in  all  cases. 

i 

i  1 

When  large  zenith  angles  are  to  be  used,  study  of  the 
temperature  profile  apd  its  time  variations  at  a  given  lo¬ 
cation  can  be  of  use  for  developing  an  N(|  profile  model  to 
,  |  fit  the  location.  The  height  parameters  for  a  theoretical 

model  based  oh  any  desired  constant  lapse  rate  are  easily 
deduced  from  the  results  of  Table  2  and  Eqs.  (11)  and  (12). 
To  match  a  giv^n  observed  integral,  the  equivalent  height 
h£j  of  the  theoretical  integral  varies  (ji  +  1 )" 1 ,  hence  the 
equivalent  height  at  03C‘  for  a  theoretical  hj  profile  of  de- 
'  gree  5  (to  fit  a  lapse  rate  of  4.8°C/km,  cf.  Fig.  1)  would 
be  (6/5)  .  40  or  48  km.  The  temperature  coefficient  of 
hj  mdst  then  chknge  by  the  same  factor.  Ray  tracing 
,  studies  may  also  be  needed, to  compare  path  curvature 
effects  for  different  temperature  profiles. 

The  larger  the  zenith  angle,  the  more  important  are 
:  i  the  effects  of  horizontal  gradients  in  the  atmosphere.  So 

far,1,  these  have  b£en  neglectfed. 

ZENITH  INTEGRAL  OF  N 

w 

i  •  « 

1  The  integral  of  Nw  (radio  value)  will  be  mentioned 
only  briefly.  For  convenience,  the  Nw  profile  can  be 
represented  by  an  expression  of  the  same  degree  as  the 
Njj  profile,  but  it  does  not  have  the  same  theoretical  ba¬ 
sis.  The  equivalent  height  hw  is  much  less  than  hj,  and 
so  far,  the  integral  of  Nw  is  predicted  from  surface  data 
witlji  an  order  of  magnitude  less  precision  than  the  dry 
integral,  i.e. ,  to  a  few  centimeters  instead  of  millimeters. 
Preliminary  work  on  Nw  was  reported  earlier  (Ref.  4), 
and  work  is  continuing  on  this  aspect  of  the  problem.  It  is 
important  for  rkdio  applications,  though  not  for  lasers. 
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Fig.  7  OBSERVED  AND  PREDICTED  VERTICAL  INTEGRAL  OF 
REFRACTIVITY  AT  WASHINGTON,  D.  C.,  JANUARY 
AND  JULY  1967,  DRY  COMPONENT 
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Fig.  8  OBSERVED  AND  PREDICTED  VERTICAL  INTEGRAL  OF 
REFRACTIVITY  AT  ALBUQUERQUE,  NEW  MEXICO, 
JANUARY  AND  JULY  1967,  DRY  COMPONENT 
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